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5_( The problem of proton-antiproton motion in the H-H system is investigated by 

c3 
^ means of the variational method. We introduce a modified nuclear interaction 

0^ through mass-scaling of the Born-Oppenheimer potential. This improved treatment 
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^_^ of the interaction includes the nondivergent part of the otherwise divergent adia- 

o batic correction and shows the correct threshold behaviour. Using this potential 

^ we calculate the vibrational energy levels with angular momentum and 1 and the 

^ corresponding nuclear wavefunctions, as well as the S-wave scattering length. We 

'""' obtain a full set of all bound states together with a large number of discretized 

^ continuum states that might be utilized in variational four-body calculations. The 

O 

l/~j results of our calculations gives an indication of resonance states in the hydrogen- 

in 

\^ antihydrogen system. 
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I. INTRODUCTION 

The interest in antihydrogen and its interaction with ordinary matter is inspired by the 
ongoing experiments on antihydrogen synthesis and trapping at CERN. The aim of these 
experiments is to use cold antihydrogen atoms for tests of the fundamental laws and sym- 
metries of Physics. 

Substantial progress has been made during the recent years in that antihydrogen atoms 
have been trapped for 1 000 s [T]. This time is sufficiently long for the H atoms to sponta- 
neously de-excite from the highly excited Rydberg states, in which they are formed, to the 
ground state that is preferred in the planned spectroscopic and ballistic experiments. 

The generic example of matter-antimatter interaction is the collision between the two 
simplest atoms of each sort, the H-H collision. In spite of its apparent simplicity the H-H 
system proves to be very challenging. It significantly differs from e.g. the H-H system by 
the presence of annihilation between particles and antiparticles and by the rearrangement: 
the impinging H and H atoms can recombine to two completely different atoms, Pn and Ps 
(see Figure [I]). The latter circumstance makes the problem very demanding both formally 
and computationally, even at the level of Coulombic description that is addressed in the 
present work. 
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FIG. 1. Two possible configurations of the system, (a) Hydrogen-antihydrogen (b) Protoniuni- 
positroniuni. 



The previous studies of the H-H system include distorted- wave approximation [2]-S] , sim- 
ple extensions thereof based on the close-coupling method |5], the optical-potential method 
[6], and the Kohn variational method [7]. All previous treatments have been based, one 
way or another, on the Born-Oppenheimer (BO) approach. The extension to the adiabatic 



treatment was not possible since it was shown that the adiabatic correction diverges [8]. 

The eigenvalues and eigenfunctions of the nuclear motion in the BO potential have so 
far been calculated only for a few of the highest exited states near the H-H dissociation 
threshold [21 EHU] , using numerical integration. In this paper we use a variational approach 
using Gaussian expansions to investigate the nuclear motion of the proton and antiproton 
in the improved BO potential. By using the variational method we obtain a set of all bound 
states (eigenvalues and corresponding eigenfunctions) in the improved leptonic potential that 
includes the non-divergent part of the adiabatic correction and shows the correct asymptotic 
threshold behaviour. 

Previous calculations indicated the presence of possible near-threshold resonance states of 
the hydrogen-antihydrogen system [2l[6l[TT]. Such states greatly influence the cross sections 
of both elastic and inelastic hydrogen-antihydrogen scattering [TTj, acting as transient states 
mediating the rearrangement to Protonium and Positronium. We therefore devote special 
care to the treatment of the near threshold states since it is a priori possible that the 
improvement of the BO potential may change the number of states below the threshold 
and/or their binding energies. The number of variationally obtained bound states below the 
threshold is double-checked using the procedure of Friedrich and Raab [121 [13] • 

The BO-approach to the H-H system works very well at the large (/im - nm) and inter- 
mediate (nm - ao) internuclear distances. However as the atoms come closer, the leptonic 
clouds start to overlap, and the H-H system gets prone to undergo a rearrangement into 
two completely different atoms, Protonium (Pn) and Positronium (Ps). The rearrangement 
is particularly probable to occur below the so called critical distance (-Re = 0.7427 [8l fll]) 
below which the proton-antiproton dipole is not able to bind the two leptons, and positron- 
ium can be released. The region around the critical distance is particularly difficult to treat. 
This is manifested by the divergence of the adiabatic correction to the Born-Oppenheimer 
interaction potential [S] and highlights the need of a full four-body treatment. 

The four-body treatment might be e.g. based on the variational method provided that 
the latter would be able to include the relevant arrangement channels and cope with their 
coupling. If the variational method utilizing an expansion in a basis is applied, the basis 
functions should preferably allow easy transformations between the various arrangement 
channels, as e.g. in the Gaussian Expansion Method (GEM) constructed by the Kamimura 
group for the four body problems [15J. Having these aspects in mind, we expand the im- 



proved BO solutions in a Gaussian basis, as to allow for the use of these solutions as a subset 
of the basis in variational 4-body calculations. We deliver full set of states that are bounded 
by the BO potential, together with a large number of discretized continuum states (all ex- 
panded in a Gaussian basis) that might be utilized in variational four-body applications. 
This might allow for a better understanding of the connection between the adiabatic and the 
full four body solutions, and provide the missing link in understanding why the adiabatic 
approach to H-H breaks down as manifested by the divergence of the adiabatic correction. 

It should be recalled that the BO energy eigenvalues and eigenfunctions are normally 
obtained in the BO potential that asymptotically converges to the H-H dissociation thresh- 
old at —1 a.u. (a value without adiabatic correction) and that was the case in all previous 
calculations. However in a 4-body calculation the H-H dissociation threshold, and all other 
thresholds, obviously takes their proper value that include the adiabatic correction. In the 
present work we correct for this defect by including the non-divergent part of the adiabatic 
correction in our modified leptonic potential through a mass-scaling of the BO potential. It 
is important to emphasize that this procedure is not equivalent to a simple shift of the thresh- 
old for the conventional BO potential (although it incidently assures the correct asymptotic 
behaviour of the modified leptonic potential). 

Atomic units have been used throughout this article. 

II. METHOD 

A. Choice of the Hamiltonian 

The 4-body Hamiltonian for hydrogen-antihydrogen system expressed in a space-fixed coor- 
dinate system reads 

o-SF ]_A ^_\ _ lA — lA 

-"4body 2mp^rp 2mp^^P 2^^c 2^^s 

+ V'(rp, Fp, re, rg) , (1) 

where rp, r^, Fc and rg are position vectors of proton, antiproton, electron and positron re- 
spectively, V describes Coulomb interactions between all particles, and rrip = 1836.15267247 
a.u. is the proton mass. Introducing body-fixed coordinates rep = re — rp, rgp = rg — rp, 
R = rp — rp and separating the center-of-mass motion, the 4-body Hamiltonian can be 



rewritten as a sum of a leptonic Hamiltonian 



iJiep = --Ar,, - -Ar,, + V{r,p, r^p, R) , (2) 

and Hamiltonian H', which in this case has the following form 

- ^Ar - TT^Ar- , (3) 

2r?ip '^cp 2mp '^cp 5 V / 

where /in = fnp/2 is a nuclear reduced mass. 

Assuming the total 4-body wavefunction as a simple product of a given leptonic function 
'^iep(i"ep,i"ep; R) , which is an eigenfunction of (pi) 

-f^lepV^lep(rep, Tep] R) = -EBo(-R)V'lep(rep, Tep] R) (4) 

and some unknown nuclear function x(R')) 0^6 obtains a nuclear Schrodinger equation within 
the adiabatic approximation 

H^XkiB.) = EkXkiB.) , (5) 

where the nuclear Hamiltonian is defined as 

H^ = -;t^Ar + ^Bo(^) + SE,^{R) . (6) 

//in 

The last term on the r.h.s of ^ is called the adiabatic correction, and is defined as an 
expectation value of the Hamiltonian H' with respect to the leptonic wave function 

SE^d{R) = (^icp(rep,rep;i?)|i7'|^iep(rep,rep;i?)) . (7) 

As it was shown by Strasburger [S], for the hydrogen-antihydrogen molecule the adia- 
batic correction diverges as the internuclear distance tends to the critical value Re ~ 0.7427 
bohr. For this reason the nuclear Hamiltonian pi) is not well defined and one cannot cal- 
culate adiabatic energy levels for this system. Of course we can apply Born-Oppenheimer 
approximation, i.e. solve Schrodinger equation with the following nuclear Hamiltonian 

H^ = -^An + E^o{R). (8) 

Z/in 

However, one can still improve the Born-Oppenheimer potential by including only the non- 
divergent part of the adiabatic correction. This can be done {vide IaI) through different 
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factorization of the 4-body Hainiltonian and can be implemented by a simple scaling of the 
Born-Oppenheimer energy. 

The last two terms on the r.h.s of (|3| do not depend on the internuclear coordinate R, 
and we can treat them as a part of the kinetic energy of leptons. This leads to a new leptonic 
Hamiltonian of the following form 

^lep = -^Arep - ^Ar,, + ^(rep, Fep, R) , (9) 

where /i = m^/{in^ + 1) is the electron-proton reduced mass. Solving the Schrodinger 
equation with Hamiltonian toj) one obtains a new leptonic energy curve ^iep(-R), which on 
top of the Born-Oppenheimer contribution includes also a part of the adiabatic correction. 
It can be shown that the new energy is related to the original Born-Oppenheimer energy by 
the following mass scaling procedure 

Ei,^{R) = fiEBoit^R) . (10) 



In can be also shown that the new leptonic potential given by (10) has the correct adiabatic 
long range asymptotic behaviour, and nonadiabatic dissociation limit. For the detailed 
derivation and discussion of this procedure see |Aj 

The scaling procedure allows us to define a new nuclear Hamiltonian with Ebo{R) + 
(5-Ead(-R) substituted by E\ep{R) 

H^ = -^An + Ei,piR), (11) 

which is well defined and can be used to calculate energy levels for the H-H molecule. This 
approach leads to the energies which should be considered as being "halfway" between the 
Born-Oppenheimer and the full adiabatic approximations. 

B. The potential energy fit 

In this work we used the Born-Oppenheimer potential for the H-H system calculated by 
Strasburger PSJ with 256 explicitly correlated Gaussian (EGG) functions. The potential 
was computed for internuclear distances R ranging from R = 0.744 bohr up to 20.0 bohrs. 
For R between 12 and 30 bohrs we have appended points computed by Strasburger with 
points obtained from asymptotic formula 

26 ^ 

Ebo{R)=E^o-J2^^ (12) 

n=6 



where E^q = — 1 hartree is the BO energy of two separated hydrogen atoms and the Van 
der Waals constants C„ were calculated by Mitroy and OvsiannikovjlTj . 
For R = Re the Born-Oppenheimer energy reaches value E^^ — 1/Rc, where Ef^ = —0.25 
hartree is the positronium ground state energy. When the internuclear distance R is smaller 
than the critical value Re, the two leptons are no longer bound in the field of a dipole 
formed by the proton and the antiproton. Therefore, for R < Re we assumed that the 
Born-Oppenheimer energy is equal to the sum of Ef^ and the Coulomb attraction between 
the nuclei 



-Sbo(-R) = -^1 



Ps 



R 



(13) 



The choice of the potential for R < Re is in the spirit of the orthodox BO- approximation, 
i.e. the potential is equal to the lowest leptonic energy at each R, which in this region 
corresponds to the positronium ground state energy. The calculated Born-Oppenheimer 
energy as well as the long and short range approximations to it are shown in Figure [2J 




FIG. 2. The Born-Oppenheimer potential (ii^Bo(-R)) ^ red, the short range approximation (13) 



green, the long range approximation ( |12| ) - magenta, and the pure leptonic potential {Ebo{R) + 
1/R) - blue. 



For our purposes we have prepared an analytical fit with 30 linear and 7 nonlinear parameters 
carefully optimized to reproduce the Born-Oppenheimer energy in all three discussed regions. 
The fitting function has the following form 



+ ELiEt=o^nkR'exp{-a^R^). 



(14) 



TABLE I. The fit parameters for Born-Oppenheimer (V^o) and mass-scaled leptonic (Mep) poten- 
tials. 





Vbo 


Mop 




Vbo 


Mop 


Aio 


-19.8582635505679 


-20.1369672678805 


A40 


-19.9771658686913 


-20.2530003255519 


All 


67.6269717956708 


44.1781330016383 


A41 


-173.8431212019852 


-131.5359580062383 


An 


-20.0575886039098 


-14.6255427563730 


A42 


39.0038819732993 


36.6202418042846 


Ai3 


1.6436298797648 


1.2854178363453 


A43 


-4.9839482825694 


2.3956839062761 


Ai4 


-0.0417677179701 


-0.0344425997356 


A44 


0.2098970274357 


0.0726049211608 


A20 


57.5162155781683 


57.6405552715681 


A50 


-22.3850547348492 


-22.6615482631176 


A21 


9.3918281802097 


-24.8917421281046 


A51 


106.8813949154074 


123.3264381665587 


A22 


3.0569545228764 


8.8916405222450 


A52 


-24.8069885152175 


-29.3520833084752 


A23 


-0.2521821480278 


-1.4152945078081 


A53 


2.2667291349955 


-4.7119434834608 


A24 


0.0965988366924 


0.1428155629355 


A54 


-0.0258580270297 


0.4403103436776 


A30 


4.7043278292101 


5.4110168738119 


Aei 


0.0000097266439 


0.0000091571743 


A31 


-16.1993647293737 


-17.2109690438145 


A62 


-0.0000006275304 


-0.0000005856539 


A32 


23.4284275566323 


24.7666406025529 


Aes 


0.0000000184890 


0.0000000171095 


A33 


-14.9456597423665 


-15.7590990169304 


A64 


-0.0000000002113 


-0.0000000001939 


A34 


3.8997649015388 


4.1982517062010 


/9 


6.1520725018366 


6.1431639772293 


01 


0.0897852714851 


0.0893701431156 


04 


0.1412060702801 


0.1171840549361 


«2 


0.2268196733512 


0.2952163755619 


"5 


0.1048123413141 


0.1108374703554 


03 


2.2437975957692 


2.2164844767807 


aa 


0.0068068098389 


0.0067006105329 



with an additional constrain 



Xl^-o 



0. 



n=l 



(15) 



The first term on the r.h.s. in (14) describes the Born-Oppenheimer dissociation hmit, the 



second one is added to ensure proper behaviour of the fit for small internuclear separations 



{vide (13)). 



We have also prepared an analytical fit for the mass-scaled leptonic potential (10). In 



this case the long range asymptotic expression for Eiep{R) has the same form as in (12), but 
with Cn constants replaced with C„ = Cn/fi"''^ and changed dissociation limit E^^. On the 



other hand the short range behaviour of Eiep{R) is derived from eqs. (13) and (10) , so the 
fitting function in this case has the following form 



+ En=i Et=o ^nkR' exp(-a„i?2) ^ 



(16) 



where E^^ = —fj, is twice the nonadiabatic ground state energy of the hydrogen atom. Also 



in this case parameters A^q are restrained by the condition 

6 
J]inO = 0. (17) 

n=l 

It should be stressed that the scahng is done for all internuclear distances, thus no disconti- 
nuity arrises in the potential (even though the critical distance changes under scaling). All 



the linear and nonlinear parameters in ( 16 ) were optimized independently of the correspond- 
ing parameters in (14). Since the values of the Born-Oppenheimer potential calculated by 
Strasburger |16j are given on a grid for chosen R: {R,E-bo{R)), we have been fitting the 
mass-scaled potential to the points obtained as follows {R/ ^, E\ep{R/ fi) = /ii?Bo(-R))- 
All parameters for (14) and (16) are given in Table III In both cases the fit errors with re- 
spect to the interaction energy are smaller than 0.1% for internuclear distances R < 8.0 



bohrs and still not larger than 0.3% for R up to 30.0 bohrs. Since functions (14) and (16) 
behave like exp{—aR^) for large R, they are not able to properly describe the asymptotic 
behaviour of the potentials for arbitrarily large R. However, we chose to use this type of 
the fitting function because it allows us to perform analytical calculations of the matrix 
elements {vide infra), which would not be possible with fits explicitly including the Van der 



Waals expansion (12). 



C. Basis functions 



To solve the nuclear Schrodinger equation with Hamiltonian defined by (pj) or (11) for a 
given angular momentum / we represent the nuclear wavefunctions as 

Xkimi^) = YiU^)MR) , (18) 

where YJ^ is a spherical harmonic function, R denotes angular coordinates of vector R, and 
(pki{R) is expressed in a basis set of analytical functions 

MR) = E ^''9^{R) = E i^2n-l9nliR) + 4Ui{R)) ■ (19) 

Basis functions g^iiR) and g^i^R) have the following form 

g'^^iiR) = N^iR' exp(-z/„i?2) cos(az/„i?2) (20) 

g:,i{R) = XiR' exp{-i^nR^) siniai^nR^) (21) 
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where a = 'n 12 and A^^;, A^^^ are normalization constants. 

The nonhnear parameters v^ defining the basis functions are chosen to be given by a 
geometrical progression 

yn = ^ (22) 



n 



r2 
with 



max 
min 



\ (n-l)/("max-l) 

r„,a. \ _ ^23) 

The parameter nmax is the number of cosine and sine oscillating Gaussian functions as defined 



in (21), used in the radial wavefunction expansion. 

Substituting ( [l8| into a nuclear Schrodinger equation one obtains the equation for the 
radial function 4>ki{R) 

^ "^^ ^ "^ +^-^!^ + V{R)-E,^MR) = 0, (24) 



2/in di?2 j^^R dR 2/in^2 

with V{R) being the Born-Oppenheimer or the mass-scaled leptonic potential. 

The linear expansion coefficients c^' defining the function (j)ki{R) are obtained variationally 
be solving the generalized eigenvalue problem 

{T + Y)cu = EuScki, (25) 

where c^; is a column vector of coefficients df , T is a kinetic energy operator matrix 

^. _/ |_^^! 1 d , /(/ + 1) , 

'' ~ ^^^' 2/i, di?2 ^^R dR + 2/ini?2 1^^/ ' i^^^ 



V is a potential energy operator matrix calculated with V{R) given by (14) or (16) 



V,, = {g,\V{R)\g^) (27) 

and S is an overlap matrix. 

For a given expansion length one may optimize the basis functions by changing the values 
of Tmin and Tmax- Howcver, we must remember that different rmin and rmax would be optimal 
for different eigenstates. Because of the shape of the potential energy, the wavefunctions 
with the lowest energies are assumed to have much smaller average radius {R) than highly 
excited wavefunctions. For this reason we need to keep the parameter rmin small enough to 
be able to describe the ground state. At the same time the parameter rmax must be large if 
we want to be able to describe the wavefunction oscillations for highly excited states. This 
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FIG. 3. Mass-scaled eigenvalues with Tj, 



0.00007 a.u. and r„ 



15 a.u. kept constant. Top: 



1/ = 26 - 28. Bottom: z^ = 29 - 32. Threshold energy E, 



lep 



-0.9994557 (red line). 



can be achieved only when the basis set expansion is long enough, which is illustrated in 
Figure |3| 



Figure |4] shows how the energies for highly excited states depend on the rmax parameter 
for fixed values of r^^^ and given expansion length. If the value of r^iax is too small the 
near threshold bound states are not described correctly. However, when rmax is too large 
one may observe oscillations of the bound states energies. We decided to use rmax = 15 for 
nmax = 60 and r^ax = 20 for rimax = 120 in our calculations. These values are large enough 
to describe all bound states but still small enough to be in the region where the oscillations 
of the bound states energies are not severe. 
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FIG. 4. The vibrational states i/ = 29 — 32 dependence on rmax, and the threshold energy (red 
line). 120 Gaussians with rmm = 0.00007 a.u. have been used. 

III. RESULTS 
A. Bound states 



In Tables [IT] and III we present the rovibrational energies of the hydrogen-antihydrogen 
molecule obtained with the Born-Oppenheimer (-E™) and the mass-scaled leptonic (-E"™) 
potentials for angular momentum / = and / = 1 respectively. With both potentials 
we have found 29 bound states for / = and 28 bound states for 1 = 1. Since for the 



small internuclear separations the potential is given by (13), which is a shifted by the Ef^ 



potential for the protonium atom, we have expected the lowest eigenvalues to be equal to 
the protonium energy plus the positronium ground state energy 



E. 



HH 



El'' + ^f ^ . 



(28) 



As it can be seen from Table [IT] this relation is well fulfilled up to the vibrational quantum 
number z/ = 20 [vide the fifth column of Table [ll]). For larger values of u the difference 
between -E™ and El'^ grows successively, which is caused by the fact that the nuclear 



wavefunction is no longer concentrated in the region where the condition ( 13 ) is fulfilled. The 



differences between the Born-Oppenheimer energy levels for I = and the protonium s states 
energies are shown in Figure [5j It is obvious that the ground state radial wavefunction for the 
protonium atom and for the H-H system are nearly identical. However, the wavef unctions of 
near threshold states of the hydrogen-antihydrogen molecule differ from protonium orbitals. 
The radial wave functions for the last two bound sates of H-H with angular momentum 
/ = are shown in Figure |6l 
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FIG. 5. The H-H oscillation levels for I = (solid blue lines) compared with protonium levels 
shifted by Ef^ (dashed red lines). The lowest shown levels correspond to quantum number u = 20. 
Black lines represent potentials: the mass-scaled BO potential (solid line), the shifted Coulomb 



potential ( 13 ) (dashed line) 
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FIG. 6. Characteristics of the bound states i/ = 28 (dashed) and i/ = 29 (solid). 240 Gaussians 



with Tn 



0.00003 a.u. and r„ 



20 a.u. have been used. 



In the two last columns in Tables [IT] and III we present the dissociation energies obtained with 
both our potentials for angular momentum / = and / = 1 respectively. The dissociation 
energy for the Born-Oppenheimer approximation is defined as 



, TTioo itHH 



whereas for the mass-scaled potential we have 






HH 



(29) 



(30) 



The difference between e^i and e^i is negligible for the low lying states. However, the differ- 
ence grows as the energies tends to the threshold energy and reaches 5% for the last bound 
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state with / = and 13% for the last bound state with / = 1. 



B. Scattering states and scattering length 



Since we are using an algebraic approximation to solve the nuclear Schrodinger equation, 
not only bound but also continuum states, however with discrete energies, are obtained (see 
Figures [T] and [s]) . Because we are using Gaussian basis functions in our description, the 
wave functions corresponding to the scattering states cannot possess the proper asymptotic 
behaviour for an arbitrarily large R. However, we were still able to provide a proper de- 
scription not only for small R, where the potential well is deep and the wave functions vary 
rapidly, but also for R large enough to see the asymptotic properties of the wavef unctions. 
To test the quality of these scattering wavefunctions we used them to estimate the value 
of the scattering length. We have used a geometrical procedure to estimate the scattering 
length, by interpreting the scattering length as an intersection of the tangent line to the 
"zero energy continuum state" R(J)e{R) with the i?-axis, where (J)e{R) is a the radial con- 
tinuum wavefunction for E ^ [18] (in our computation the lowest discretized continuum 
state, see Figure |9|). 
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FIG. 7. The first few continuum states z^ = 30 (dashed), z/ = 31 (solid) and z/ = 32 (dotted). 240 
Gaussians with r^i^i = 0.00003 a.u. and r^ax = 20 a.u. have been used. 
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FIG. 8. The continuum states z^ = 33 (solid), z^ = 34 (dashed) and v = 35 (dotted). 240 Gaussians 
with Tinin = 0.00003 a.u. and rmax = 20 a.u. have been used. 
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FIG. 9. The last bound state i/ = 29 (solid), the first continuum state i/ = 30 (dashed), and the 
tangent line to the first continuum state in the large R asymptotic region (red). 240 Gaussians 
with rjnin = 0.00003 a.u. and rmax = 20 a.u. have been used, and the estimated scattering length 
is Rsc ~ 7.6. 

The estimated scattering length obtained with this procedure is a = 7.6±0.4 bohr, where the 
uncertainty conies from the sensitivity of the obtained results to the Gaussian parameters 
Tinax and nmax- This result can be compared with the value a = 7.7 bohr, which was obtained 
by numerical integration of the Schrodinger equation [llj. This proves that our calculations 
have obtained a reasonably good description of the scattering states wavefunctions. 

C. Relation to the WKB approach 



The semi-classical approximation has been proved to give accurate predictions of the number 
of bound states as well as the value of the scattering length for diatomic molecules [T2| 
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We have decided to use the WKB approach to confirm the results we have obtained with 
the algebraic approximation. 

When the potential energy has an attractive tail which decays faster then l/R"^, the 
system has a finite number of bound states where the dissociation energies e^ and quantum 
numbers i/ = 1, 2, 3, . . . obey the following quantization rule 

z/th - z/ = F(e,) , (31) 

where z/^h is, in general non-integer, threshold quantum number, and F is called the quanti- 
zation function. As it was shown by Friedrich and Raabp^ [T3] for a potential well with a 
homogeneous tail of order (—6), the F function is given in the following form 



F{e. 



X 



27r[l+(«;,/36)4] "^ 1+(k,/36)4 

1 , D , r(|)K^2/3 

7^- 



+ n^l., a„^2/3 + 



(32) 



8 ^ 27r(K,/36)2/3 ^ 40Fr(|) 

where the constants 6, d, B and D are defined in Table I in fT2] , k^, is defined as n^ = y/2jAe^, 
where A^ is a nuclear reduced mass and Pq is given in terms of the strength of the van der 
Waals interaction Ce as 

13, = {2CeMY^^ . (33) 

The Pq defines the potential range and is a typical scale for the quantum mechanical wave- 
lengths and penetration depths. In our case Pq = 10.4521 for the Born-Oppenheimer po- 
tential and Pq = 10.4592 for the mass-scaled potential. We have calculated values of the 
F(e^) function for the dissociation energies listed in Table [11} Substituting these results into 



formula (31) we obtained the threshold quantum numbers z/th. The values of z/th computed 



from subsequent dissociation energies e^, are shown in Table IV and Figure 10 The predicted 
number of the bound states is the largest integer less than z/th- For lower quantum numbers 
z/ the predicted number of bound states is rapidly changing, but for z/ > 22 it stabilizes at a 
value equal 29, with one exception for z/ = 28 where the z^th value is slightly below 29. The 
obtained results prove that we have succeeded to find all the bound states for H-H molecule 
in its leptonic ground state and angular momentum zero. 



Using the quantization function (32) and the relation 



a H , ^, ,, (34) 

tan(7rF(e,)) ^ ^ 
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we were able to obtain an independent estimation of the scattering length. In the above 
formula a = b = 0.4779888/36 (see Table I in [I2])- Substituting dissociation energies of 
the last bound state (zero angular momentum) calculated with the BO potential and with 
the mass-scaled BO potential, one gets the scattering length equals 7.6 bohr and 7.5 bohr 
respectively. These values are in full agreement with the results obtained from the alge- 
braic approximation, which proves that we succeeded to correctly describe the asymptotic 
behaviour of the scattering wavefunctions in our calculations. 



TABLE IV. The threshold quantum number obtained for the BO potential (fth) and the mass- 
scaled BO potential (t'th)- 



I'th 



i^th 



20 


28.75 


28.76 


21 


28.92 


28.93 


22 


29.07 


29.08 


23 


29.20 


29.21 


24 


29.30 


29.31 


25 


29.34 


29.35 


26 


29.29 


29.29 


27 


29.07 


29.08 


28 


28.94 


28.94 


29 


29.34 


29.35 




FIG. 10. The Quantization function for the BO potential for states with quantum number v 
6-29. 
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IV. CONCLUSIONS 

By applying the algebraic approximation we have succeeded in obtaining the vibrational 
spectra of the leptonic ground state of the hydrogen-antihydrogen molecule in its two lowest 
rotational states. Since there is no possibility to include the adiabatic correction for this 
system, our calculations were done at the Born-Oppenheimer level. However, we were able 
to append the Born-Oppenheimer potential with the nondivergent part of the otherwise 
divergent adiabatic correction. This contribution is responsible for the change in the long 
range behaviour of the interaction potential, it mimics the total adiabatic correction for large 
internuclear distances R and converges to the correct dissociation threshold. The discussed 
part of the adiabatic correction was included via the mass-scaling procedure, which does not 
require any additional calculations besides the Born-Oppenheimer one. The mass-scaling 
procedure leads to the proper nonadiabatic dissociation threshold of the leptonic potential. 
As it was shown the near threshold states are significantly affected by this procedure, since 
the changes in the dissociation energies reaches a few percents (as much as 13% for the most 
loosely bound state of the / = 1 symmetry). 

Using the algebraic approximation we were able to obtain wavefunctions not only for all 
bound states, but for continuum states as well. The energies of the calculated scattering 
states are of course discretized. However, the corresponding wavefunctions reproduce the 
correct asymptotic behaviour up to a certain internuclear distance R* that is larger than the 
range of the spacially extended near-threshold states. This allowed us to calculate the scat- 
tering length for the hydrogen-antihydrogen collisions. The value obtained in this calculation 
is in very good agreement with the results of high precision numerical computations and 
with the results of semi-classical approach, which confirms good quality of the discretized 
continuum states. 

The existence of the bound states for the H-H system in the Born-Oppenheimer approx- 
imation may indicate the existence of quasi-stable states (resonances) for this system in the 
nonadiabatic description. To confirm or falsify these predictions one has to perform 4-body 
calculations for this system. The wave functions obtained in our calculations can be used 
as a well adapted basis set for the internuclear degree of freedom in the H-H channel in 
the fully nonadiabatic (4-body) calculations for H-H. We suggest that the nuclear wave 
functions obtained in the present work will provide a better description of the H-H channel 
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in the 4-body calculations as compared to the alternative description using protonium basis 
set. Application of the here obtained nuclear wave functions as a part of the basis set in 
4-body calculations may not only be more efficient but also allow for better understanding 
of the relation between the B-0 and fully non-adiabatic methods in the demanding case of 
H-H. Using a mass-scaled BO basis will not cure the nonadiabaticity, but it might facilitate 
its investigation. The fact that the BO functions are part of the full basis helps to iden- 
tify the adiabatic components in the nonadiabatic solutions. In the 4-body calculations the 
mass-scaled BO basis must be augmented by other dedicated subspaces for description of 
other Jacobi fragments and coordinates. 
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Appendix A: Mass-Scaling procedure of the Born-Oppenheimer potential 

The adiabatic correction to the BO potential is defined as an expectation value of the 
nuclear part of the Hamiltonian with respect to the leptonic wave function (i.e. the eigen- 
function of the leptonic Hamiltonian). As it was shown by Strasburger [H] in the case of 
hydrogen-antihydrogen interaction the adiabatic correction is not well defined for certain in- 
ternuclear separations, and the standard procedure of improving the BO potential cannot be 
applied. For this reason we decided to apply another approach, which enables us to include 
at least the part of the adiabatic correction which is responsible for long range behaviour 
of the interaction potential. This is done following a different definition of the leptonic 
Hamiltonian ^U\ 120] and results in a simple scaling of the Born-Oppenheimer potential. 

The total Hamiltonian for the H-H molecule expressed in a space-fixed coordinate system 
reads 

Htot = Tiep + T^nc + V , (Al) 
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where 

Tiep = -^Are-^Arg, (A2) 

1 1 
T„uc = Ar Ar- , (A3) 

and V describes Coulomb interactions between all the particles in the system. 

The first step in most rigorous treatments is the separation of the COM motion. The 
COM coordinates are defined as 

X = -— (rc + rg + mprp + mprp) , (A4) 

where M = 2mp + 2 is the total mass of the system. With this separation made, the total 
Hamiltonian can be rewritten as a sum of the COM kinetic operator and the Hamiltonian 
for the internal motion of leptons and nuclei 

i/tot = -^Ax + i/rcl. (A5) 

The total wave function separates as 

^tot(rc, rg, Fp, Fp) = ^rc^r^;., • • • )x(X) , (A6) 

where r^,y are body-fixed coordinates. Following the Born-Oppenheimer approximation, the 
Hamiltonian describing the relative motion is divided into Hi^p and H', where the leptonic 
Hamiltonian Hi^p does not depend on the nuclear masses 

ifrel = -f^lep + H . (A7) 

If the internal coordinates are chosen to be 

I"ep I'c I"p ) I'ep I"e I"p ; -TV, Fp Tp , V-^^/ 

the leptonic part reads 

Hlep = -^ Ar,, - ^ Ar,, + ^(rep, Fgp, R) (A9) 

and the H' operator itself can be divided into three contributions 

Hi. = -^An, (AlO) 



rrip 



i^(',) = -^Vk (Vr., - Vr J , (All) 



mp 



H'f , = ^Ar ^Ar- . (A12) 

(^) 2mn """ 2mn "'" ^ ' 
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Part (a) describes the relative motion of the nuclei, (b) couples the nuclear and leptonic 
motion, and (c) can be considered as a correction to the leptonic kinetic energy. In fact Hi-. 
has the same form as the kinetic energy of leptons, but it is divided by the proton mass. 



The factorization of Hrei ( A7) given by (A9) and (A12) is not unique, and we can choose 
to include the H',^-. contribution in the leptonic Hamiltonian 



-f^lep — Hicp + i/(c) 



Ar 

2/i '" 



2/i 



Ar,p + V{r,p, rep, R) . 



(A13) 



This leads to the operator of the same form as (A9), but with the electron mass replaced by 



the electron-proton reduced mass, /i = nip/ {nip + 1). 



One can relate the old leptonic Hamiltonian (A9) to the new one (A 13) by applying the 
following coordinates transformation 



•■ep 



Sep /^I'ep 



Tep ' Sep fJj'^ep 



R 



S = /iR. 



(A14) 



With these definitions we can write 



-^lepl^Sep, Sep, oj 



r)— 'Sep ' r)^Sop + I' (Sep, Sep, oj 



-A, 
i 

1 



2/^2 



^Ar,p + -\^(rep,rep,R), 



(A15) 



where we made use of the fact that the potential energy operator l^(rep, rgp, R) is, in all 
coordinates, a homogeneous function of order —1. This leads to the following relation 



-f^lep(l"ep, rgp, R) — P-HiepifJ^rep, ^Yep, A'-R) 



(A16) 



It is therefore clear that the eigenvalues and eigenvectors of the new leptonic Hamiltonian 

(A17) 



(A13) 



-f^lepV^leplfep, ^ep; R) = -Elep(-R)V'lep(l'ep, ^ep; R) 



and the eigenvalues and eigenvectors of the original leptonic Hamiltonian (A9) 



Hlepl/Jlepirep, fep; R) = -£^Bo(-R)V'lep(l"ep, ^ep; R) 



can be related by the following scaling transformations 



(A18) 



^lep(i?) = ^E^oi^R) 



(A19) 
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and 



V'lep(rep, rgp; R) = /i^^lep(yUrep, flT^p] flR) 



(A20) 



where /i^ on the r.h.s of the above equation is a normahzation factor. 



It is obvious that the mass-transformed leptonic energy (A19) consist of the Born 



Oppenheimer energy Ebo{R) and some non-Born-Oppenheimer contributions. To discuss 
this in a systematic way, let us express E\ep as an expectation value of the H\ep operator 
with the wavefunction ipi^p 



Eiep(i?) = (V'lep(rep, ^ep] i?) |i?lep|^lep(rep, ^ep] R)) 



Using (A20) we get 



2cW/ 



E,,p{R) = fi'TifiR) + fiVifiR) + fi'S'^^ifiR) 



(A21) 



(A22) 



where T and V denotes the expectation values of the kinetic and potential energy operators 
respectively 



and 



1 1 

T{R) = (^iep(rep,rep;i?)| - -Ar,p - -Ar,p|^icp(rcp,rcp;i?)) , 

V{R) = (V'iep(rcp,rep;/?)|V(rep,rep,R)|V'iep(rep,rep;i?)) 



^id(^) = (V^lep(rep,rep;/?)|i^(c)|V^lep(rep,rep;i?)) 



(A23) 
(A24) 

(A25) 



is a part of the total adiabatic correction £ad = (V^l-f^'l"^)- By expanding (A22) in powers of 
m~^ one obtains the following expression 

Eiep(i?)= Ej,oiR)+£i^iR) 

- '2T{R) + V{R) + R^-^^] i- + 0{^) . (A26) 



On the other hand, applying eqs. (A16), (A19), (A20) and the Hellmann-Feynmann theorem 
we can write 



(^lep (/iFep , fiTep] fJ,R)\— i^lep (/iFep , /Xrgp , /iR) | Iplep (/iF ep , /irgp ; /ii?) ) 



=1 dfi 



EBoifiR) 



(A27) 



Once again using the fact that l^(rep, rgp, R) is a homogeneous function of order —1, one 
obtains the following form of the virial theorem [2T1 [22] 

2r(i?) + V(i?) + r'^^^^ = , (A28) 

dR 



23 

which proves that the expression in the square brackets in (A26) equals zero. Finally £'iep(i?) 



reads 

Eiep(i?) = Ebo{R) + £if{R) + 0{m;^) . (A29) 

So the mass-scaled leptonic energy E{R) consists of the Born-Oppenheimer energy E{R), 
a part of the adiabatic correction £l^{K) and some other higher order corrections, which 
should be classified as nonadiabatic contributions. These nonadiabatic corrections can be 
written in the following explicit form 

-Elep(-R) - -Ebo(-R) - ^ad (R) = 

= Y.n=2 (^j ELo (fl) (;^dk^Bo(^) • (A30) 

The adiabatic effects included by the mass-scaling procedure come from the monomer 
part of the H' operator. The two other parts (given as expectation values of H',-, and H',^-, 
operators with the leptonic wave function) are not included in the E(R), however, the E^^ 
correction is important to correctly reproduce the asymptotical behaviour of the adiabatic 
potential. For large internuclear separations R the interaction part of the Born-Oppenheimer 
energy decays as 

Ebo{R) - ^Bo(oo) = -C^R-^ - CsR-^ - CioR-'' + • • • (A31) 



Due to the scaling relation (A19) the interaction part of the mass-scaled leptonic energy for 
large R behaves as 

Eiep(i?) - ^iep(oo) = -C,R-' - CsR-' - C,oR-'' + • • • 

= - ^i?-6 - ^i?-8 - ^i?-io + ■ • • (A32) 

Expanding Cq in powers of m~^ one gets 

C, = C, + —C, + 0{m-^) = C, + 5Cf + 0(m-2) , (A33) 

r/ip ^ ^ 

where dC^'^ = SCg/mp is the adiabatic correction to the Cq constant [23j. 

It should also be stressed that after the scaling procedure the leptonic energy curve has 
a proper nonadiabatic dissociation limit 

lim Eiep(i?) = -/i, (A34) 
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i.e. twice the nonadiabatic energy of the ground state of the hydrogen atom. This is a 
simple consequence of using Hi^p instead of Hiep since the former contains the nonadiabatic 
Hamiltonian for the monomers. 

In Table |V] we compared the values of the adiabatic correction £^ad calculated by Stras- 
burger [8] with the correction obtained from the scaling procedure 6E\ep{R) which is the 
difference between the mass-scaled and the Born-Oppenheimer energy 

SEi,p{R) = Eiep(i?) - Ebo(^) (A35) 

and A(_R) is defined as 

A{R) = SEi,p{R) - S^{R) . (A36) 

For large internuclear separations A(i?) tends to the difference between the nonadiabatic 
and the adiabatic dissociation threshold (i.e. —1 + l/mp) 

lim A(R) = -/i + 1 = . (A37) 

i?-s>oo rrip nip + 1 nip 

This difference comes from the nonadiabatic contributions 0{nip^) which are included in 



(A29). 



As it can be seen from Table |V] the scaling procedure reproduces over 95% of the total 
adiabatic correction for R larger than 2.0 bohrs. For smaller internuclear distances the in- 
teraction part of the adiabatic correction becomes larger and the difference grows. However, 
in this region the adiabatic approximation starts to fail, and as R gets close to the critical 
distance R^ the adiabatic correction diverges. Since the S^^ correction included in our 
scaled energy does not diverge for R ^ Re, one can conclude that the other contributions 
are responsible for this divergent behaviour of the adiabatic correction near the critical dis- 
tance. This conclusion seems to be in a good agreement with the reasoning presented by 
Strasburger. It is shown in p^j that with a rough approximation the expectation value of 
the Vr operator with the leptonic wavefunction is proportional to the mean value of the 
distance between leptons and nuclei. As R tends to -Re the leptonic wavefunction becomes 
more and more diffused due to the fact that leptons become more weakly bound. In the 
limit when R reaches the critical value, the leptons are no longer bound to the nuclei and the 
mean value of the distance between them is infinite. Since the H',^. operator does not contain 
a differentiation over the internuclear distance R, the £J^'^ part of the adiabatic correction 
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is not affected by the mechanism described above and is well defined for any R. The other 
two contributions, i.e., (Hi-.) and {Hi^.), are supposed to diverge as R tends to Re- 
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TABLE II. The energy levels of the H-H molecule for angular momentum I = 0: E^q - obtained 
with the Born-Oppenheimer potential (second column); -E™ - obtained with the mass-scaled 
potential (third column). E^^ - energy levels for the protonium atom (fourth column); 6u = 
El"^ — E^Q - Difference between protonium and H-H energies (fifth column); e^Q - dissociation 
energies obtained with the Born-Oppenheimer potential (sixth column); e^o dissociation energies 
obtained with the mass-scaled potential (seventh column). 



u 


^.T 


pHH 


£P" 


s. 


<^vO 


fj/O 


1 


-459.28810584 


-459.28797088 


-459.03816812 


0.250 


458.28810584 


458.28851520 


2 


-115.00950131 


-115.00936938 


-114.75954203 


0.250 


114.00950131 


114.00991370 


3 


-51.25422639 


-51.25409833 


-51.00424090 


0.250 


50.25422639 


50.25464265 


4 


-28.93998669 


-28.93986212 


-28.68988551 


0.250 


27.93998669 


27.94040644 


5 


-18.61185049 


-18.61172801 


-18.36152672 


0.250 


17.61185049 


17.61227233 


6 


-13.00167067 


-13.00154836 


-12.75106023 


0.251 


12.00167067 


12.00209268 


7 


-9.61898924 


-9.61886541 


-9.36812588 


0.251 


8.61898924 


8.61940973 


8 


-7.42343727 


-7.42331085 


-7.17247138 


0.251 


6.42343727 


6.42385517 


9 


-5.91797886 


-5.91784964 


-5.66713788 


0.251 


4.91797886 


4.91839396 


10 


-4.84088720 


-4.84075555 


-4.59038168 


0.251 


3.84088720 


3.84129987 


11 


-4.04378975 


-4.04365618 


-3.79370387 


0.250 


3.04378975 


3.04420050 


12 


-3.43753603 


-3.43740095 


-3.18776506 


0.250 


2.43753603 


2.43794527 


13 


-2.96591308 


-2.96577681 


-2.71620218 


0.250 


1.96591308 


1.96632113 


14 


-2.59194272 


-2.59180584 


-2.34203147 


0.250 


1.59194272 


1.59235016 


15 


-2.29036921 


-2.29023270 


-2.04016964 


0.250 


1.29036921 


1.29077702 


16 


-2.04344058 


-2.04330543 


-1.79311784 


0.250 


1.04344058 


1.04384975 


17 


-1.83852260 


-1.83838880 


-1.58836736 


0.250 


0.83852260 


0.83893312 


18 


-1.66669147 


-1.66655687 


-1.41678447 


0.250 


0.66669147 


0.66710119 


19 


-1.52177984 


-1.52163968 


-1.27157387 


0.250 


0.52177984 


0.52218400 


20 


-1.39961222 


-1.39945939 


-1.14759542 


0.252 


0.39961222 


0.40000371 


21 


-1.29734882 


-1.29717481 


-1.04090288 


0.256 


0.29734882 


0.29771913 


22 


-1.21293458 


-1.21273047 


-0.94842597 


0.265 


0.21293458 


0.21327479 


23 


-1.14468025 


1.14443786 


-0.86774701 


0.277 


0.14468025 


0.14498218 


24 


-1.09104749 


-1.09075955 


-0.79694126 


0.294 


0.09104749 


0.09130387 


25 


-1.05068830 


-1.05034841 


-0.73446107 


0.316 


0.05068830 


0.05089273 


26 


-1.02258176 


-1.02218304 


-0.67905054 


0.344 


0.02258176 


0.02272736 


27 


-1.00600201 


-1.00553517 


-0.62968199 


0.376 


0.00600201 


0.00607949 


28 


-1.00065727 


-1.00012002 


-0.58550787 


0.415 


0.00065727 


0.00066434 


29 


-1.00004946 


-0.99950777 


-0.54582422 


0.454 


0.00004946 


0.00005209 


30 


-0.99999650 


-0.99945226 


-0.51004241 


0.490 


-0.00000350 


-0.00000342 
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TABLE III. The H-H energy levels and dissociation energies for angular momentum / = 1: the 
energy levels E^^ and -E™ obtained with the Born-Oppenheimer and mass-scaled leptonic poten- 
tials respectively; the dissociation energies e^i and e^i obtained with the Born-Oppenheimer and 
mass-scaled leptonic potentials respectively. 



u 


pHH 


^.T 


tul 


El/l 


1. 


-115.00950985 


-115.00937728 


114.00950985 


114.00992160 


2. 


-51.25422287 


-51.25409434 


50.25422287 


50.25463866 


3. 


-28.93997500 


-28.93985016 


27.93997500 


27.94039448 


4. 


-18.61183446 


-18.61171187 


17.61183446 


17.61225619 


5. 


-13.00165510 


-13.00153280 


12.00165510 


12.00207712 


6. 


-9.61897877 


-9.61885498 


8.61897877 


8.61939931 


7. 


-7.42343471 


-7.42330837 


6.42343471 


6.42385269 


8. 


-5.91798413 


-5.91785500 


4.91798413 


4.91839932 


9. 


-4.84089719 


-4.84076562 


3.84089719 


3.84130994 


10. 


-4.04379959 


-4.04366610 


3.04379959 


3.04421042 


11. 


-3.43754139 


-3.43740635 


2.43754139 


2.43795067 


12. 


-2.96591238 


-2.96577614 


1.96591238 


1.96632046 


13. 


-2.59193806 


-2.59180118 


1.59193806 


1.59234550 


14. 


-2.29036510 


-2.29022859 


1.29036510 


1.29077291 


15. 


-2.04344083 


-2.04330568 


1.04344083 


1.04385000 


16. 


-1.83852684 


-1.83839303 


0.83852684 


0.83893735 


17. 


-1.66669349 


-1.66655892 


0.66669349 


0.66710324 


18. 


-1.52176867 


-1.52162862 


0.52176867 


0.52217294 


19. 


-1.39957504 


-1.39942243 


0.39957504 


0.39996675 


20. 


-1.29727458 


-1.29710087 


0.29727458 


0.29764520 


21. 


-1.21281710 


-1.21261340 


0.21281710 


0.21315772 


22. 


-1.14452007 


-1.14427809 


0.14452007 


0.14482241 


23. 


-1.09085090 


-1.09056341 


0.09085090 


0.09110773 


24. 


-1.05046575 


-1.05012623 


0.05046575 


0.05067055 


25. 


-1.02234924 


-1.02195075 


0.02234924 


0.02249507 


26. 


-1.00579048 


-1.00532349 


0.00579048 


0.00586781 


27. 


-1.00056446 


-1.00002650 


0.00056446 


0.00057082 


28. 


-1.00001426 


-0.99947180 


0.00001426 


0.00001612 


29. 


-0.99999527 


-0.99945098 


-0.00000473 


-0.00000470 
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TABLE V. Comparison of the non-Born-Oppenheimer correction 
procedure {5Eiep{R)) and the adiabatic correction calculated by 
milihartrees. 



obtained by the mass-scaling 
Strasburger [8j {£ad{R)) ~ in 



R 


SEi^piR) 


£.d{R) 


A(R) 


R 


5Ei,p{R) 


£.diR) 


A{R) 


0.744 


0.1441275 


5934.3 


-5934.2 


2.7 


0.5395071 


0.5328086 


0.0066985 


0.746 


0.1445682 


6513.1 


-6513.0 


2.8 


0.5410885 


0.5338338 


0.0072547 


0.748 


0.1450162 


2067.8 


-2067.7 


2.9 


0.5422054 


0.5347878 


0.0074176 


0.75 


0.1454715 


1469.23 


-1469.08 


3.0 


0.5429322 


0.5356682 


0.0072640 


0.8 


0.1591394 


40.385 


-40.226 


3.1 


0.5433773 


0.5364731 


0.0069042 


0.85 


0.1767098 


11.8294 


-11.6527 


3.2 


0.5436549 


0.5372086 


0.0064463 


0.9 


0.1972185 


5.656181 


-5.458963 


3.3 


0.5438564 


0.5378795 


0.0059769 


0.95 


0.2195822 


3.36107 


-3.14149 


3.4 


0.5440330 


0.5384894 


0.0055436 


1.0 


0.2427954 


2.27030 


-2.02750 


3.5 


0.5441933 


0.5390447 


0.0051486 


1.1 


0.2887498 


1.310081 


-1.021331 


3.6 


0.5443153 


0.5395488 


0.0047665 


1.2 


0.3311600 


0.9238342 


-0.5926742 


3.8 


0.5443107 


0.5404270 


0.0038837 


1.3 


0.3685480 


0.7397641 


-0.3712161 


4.0 


0.5438799 


0.5411558 


0.0027241 


1.4 


0.4005566 


0.6437032 


-0.2431466 


4.2 


0.5432116 


0.5417607 


0.0014509 


1.5 


0.4273892 


0.5908774 


-0.1634882 


4.5 


0.5426022 


0.5424800 


0.0001222 


1.6 


0.4496878 


0.5611774 


-0.1114896 


5.0 


0.5436598 


0.5433068 


0.0003530 


1.7 


0.4683043 


0.5444602 


-0.0761559 


5.5 


0.5444552 


0.5438155 


0.0006397 


1.8 


0.4839570 


0.5352240 


-0.0512670 


6.0 


0.5438428 


0.5441254 


-0.0002826 


1.9 


0.4970463 


0.5304100 


-0.0333637 


6.5 


0.5436756 


0.5443113 


-0.0006357 


2.0 


0.5077419 


0.5282176 


-0.0204757 


7.0 


0.5440693 


0.5444247 


-0.0003554 


2.1 


0.5162030 


0.5275680 


-0.0113650 


8.0 


0.5442068 


0.5445352 


-0.0003284 


2.2 


0.5227266 


0.5278091 


-0.0050825 


9.0 


0.5442654 


0.5445789 


-0.0003135 


2.3 


0.5277358 


0.5285403 


-0.0008045 


10.0 


0.5442959 


0.5445979 


-0.0003020 


2.4 


0.5316630 


0.5295251 


0.0021379 


12.0 


0.5443144 


0.5446109 


-0.0002965 


2.5 


0.5348346 


0.5306166 


0.0042180 


15.0 


0.5443191 


0.5446156 


-0.0002965 


2.6 


0.5374283 


0.5317269 


0.0057014 


20.0 


0.5443203 


0.5446168 


-0.0002965 



